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Abstract 

The production of invisible particles plays great importance in high energy physics. Large part of interesting 
electroweak processes include production of neutrinos, while many new physics scenarios predict the existence of 
similarly weakly-interacting particles. In events with associated production of invisible particles and hadronic jets, 
the measurement of the imbalance in transverse momentum of the final state particles is the major leverage to reject 
the otherwise dominant source of backgrounds in hadron colliders, i.e. the generic production of many jets by QCD 
interactions. Here we discuss a novel technique which utilizes the information derived from the spectrometer, even- 
tually coupled with the more straightforward calorimeter information, to infer the passage of invisible particles. We 
check the validity of this technique in data and Monte Carlo simulations in a broad range of topologies, starting from 
the simplest, with two jets in the final state, to the ones with very large jet multiplicities. We also suggest a new 
way, based on the same approach, to measure the yields and model the kinematics of the QCD multijet background in 
invisible particles plus jets signatures. The results are derived using data collected with the CDF II detector; we argue 
that the application to LHC experiments is straightforward. 
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1. Introduction 

The associated production of invisible particles and 
hadronic jets is distinctive of a large number of the top- 
ical scenarios currently investigated at hadron colliders . 
The prototypical particle that goes undetected in col- 
lider experiments is the neutrino. Neutrinos, which are 
produced in large part of the electroweak processes, in- 
teract only through the weak force, and their passage is 
inferred indirectly from the imbalance of the total trans- 
verse energy of the final state particles. Final states in- 
cluding neutrinos and quarks are of large interest in a 
variety of measurement or searches; a notable example 
is the search for the associated production of an Higgs 
boson and a Z boson, in events where the Higgs boson 
decays to a bb pair, and the Z boson decays to two neu- 
trinos til|2[]. 

Many physics scenarios beyond the Standard Model 
predict the existence of particles with such a low in- 
teraction cross section to result in an imbalance in 
total transverse energy. Among these, the quest for 
dark matter is a crucial topic in high energy physics 



field. Searches at colliders can take advantage of jets 
from initial state radiation to explore direct dark mat- 
ter production [3]. Recent experimental work suggest 
that these kind of searches can provide complementary 
information with respect to direct experiments on the 
mass and cross section of the hypothetical dark mat- 
ter particles fl, ^ . Another example is Supersymme- 
try, which provides a whole new spectrum of particles, 
the lightest of which have to be stable, neutral and non 
strong-interacting. Events with production of squarks or 
gluinos are expected to include strong-interacting SM 
particles in the final state, to conserve the color charge, 
together with lightest supersymmetric particles, result- 
ing in missing transverse momentum plus jets signa- 
tures. 

In all instances, the presence of undetected particles 
is signaled in collider detectors by an imbalance in the 
transverse momentum of the final products, which is 
generally measured with the use of calorimeters and is 
traditionally referred to, in literature, as "missing trans- 
verse energy" (¥t)- The presence of large missing trans- 
verse energy is a distinctive feature of all the signatures 
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we have cited, allowing to significantly reduce the con- 
tribution of the generic production of strong-interacting 
particles through the more mundane quantum chromo- 
dynamics, (i.e. QCD background). Still, the presence 
of QCD is unavoidable due to calorimeter resolution ef- 
fects. 

The challenging and hostile environment for invisible 
particle detection at hadron colliders thus calls for the 
need of innovative detection procedures. In this work 
we analyze the performance of a novel invisible particle 
detection technique, which is based on track informa- 
tion for evaluation of undetected particles' momentum. 

New techniques for this purpose are needed in par- 
ticular for kinematics in which the invisible particles in 
the final state carry low transverse momentum, because 
this imply low ¥t- For example, one of the last super- 
symmetric scenario not ruled out by LHC data is the 
one where the next-to-lightest supersymmetric particle 
is one of the 3rd generation squarks, the stop (J) and 
sbottom (b), the lightest being the neutralinos |6]. After 
their pair production, t (b) would likely decay to a top 
(bottom) quark and a neutralino x mat g° es undetected. 
This signature is particularly difficult to investigate if 
the difference between Mt(M~ b ) and M t + M x (M b + M x ) 
is small, because in this case the decay products have 
low momentum. The fact that experimentalists tend to 
use very stringent cuts on the event missing transverse 
energy, in order to reduce the QCD background, drasti- 
cally reduce sensitivity in this region of the SUSY pa- 
rameters space (see [7] as an example). 

Dark Matter searches face the same problem if the 
mass of the Dark Matter particle is >100 GeV/c 2 be- 
cause in this case large fraction of the collisional en- 
ergy becomes mass. Together with the lower production 
cross section, this makes collider experiments less sen- 
sitive than direct detection experiments for Dark Matter 
masses higher than few tens of GeV/c 2 , even in the case 
in which the direct detection bounds are limited by a 
spin-dependent cross section hypothesis 10, 0]. The use 
of QCD rejection techniques like the one we describe in 
the following could help extend experimental limits to 
these unexplored regions, since collider experiments are 
much less limited by the spin-dependent hypothesis. 

In all these cases, since the expected ¥t is relatively 
low, square cuts on at trigger level would not succeed 
to efficiently reject QCD background, maintaining at the 
same time high efficiency on the signal. Therefore, new 
techniques to identify, reject and eventually model the 
QCD background are needed to extend limits in these 
regions of the phase space. 



2. The CDF II detector 

The CDF II detector has been extensively described 
elsewhere in the literature fl. The detector is cylindri- 
cally symmetric around the pp beam line; a magnetic 
spectrometer consisting of tracking devices inside a 3-m 
diameter, 5-m long superconducting solenoidal magnet 
with an axial magnetic field of 1 .4 T measures the mo- 
menta and trajectories of charged particles. A five layer 
double-sided silicon microstrip detector (SVX) covers 
the region between 2.5 to 11 cm from the beam axis. 
Three separate SVX barrel modules along the beam line 
cover a length of 96 cm, approximately 90% of the lu- 
minous beam interaction region. Three of the five lay- 
ers combine an r-(p measurement on one side and a 90° 
stereo measurement on the other, and the remaining two 
layers combine an r-0 measurement with small angle 
stereo at ±1.2°. The typical silicon hit resolution is 11 
yum. Additional Intermediate Silicon Layers (ISL) at 
radii between 19 and 30 cm from the beam line in the 
central region link tracks in the SVX to hits in the Cen- 
tral Outer Tracker (COT), a 3.1-m long open-cell drift 
chamber [9] which occupies the radial range 40-137 cm. 
Eight superlayers of drift cells with 12 sense wires each, 
arranged in an alternating axial and ±2° pattern, provide 
up to 96 measurements of the track position. The po- 
sition resolution of a single drift time measurement is 
about 140 jim. Full radial coverage extends up to I77I < 
1 for the COT and up to \rj\ < 2 for the silicon detectors. 

Outside the tracking systems and the solenoid, 
segmented calorimeters with projective geometry are 
used to reconstruct electromagnetic (EM) showers and 
jets. The calorimeter is segmented radially into lead- 
scintillator electromagnetic and iron- scintillator 
hadronic ifllll compartments, both segmented into tow- 
ers, each covering a small range of pseudorapidity and 
azimuth, and in full cover the entire In range in az- 
imuth and the pseudorapidity regions of \tj\ <1.1 (cen- 
tral) and 1.1< I77I < 3.6 (plug). The transverse en- 
ergy E T = EsinO, where the polar angle 6 is calcu- 
lated using the measured z position of the event ver- 
tex, is measured in each calorimeter tower. Proportional 
and scintillating strip detectors measure the transverse 
profile of EM showers at a depth corresponding to the 
shower maximum. Electrons are identified in the cen- 
tral EM calorimeter as isolated, mostly electromagnetic 
clusters that also match with a track in the pseudora- 
pidity range \r]\ < 1.1. The electron transverse energy is 
reconstructed from the electromagnetic cluster with pre- 
cision 5(E T )/E T = n.5%HE T (GeV)@ 2%, where the 
© symbol denotes addition in quadrature. Jets are iden- 
tified as electromagnetic and hadronic calorimeter clus- 
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ters using the jetclu algorithm Ill2ll with a cone size of 
0.4, and their energies are corrected for the calorimeter 
non-linearity, losses in the gaps betwen towers, multiple 
primary interactions, underlying event and out-of-cone 
losses II 1 311 - 

Drift chambers located outside the central hadronic 
calorimeters and behind a 60 cm thick iron shield de- 
tect muons with \rj\ < 0.6, while muons in the re- 
gion between 0.6 < \rj\ < 1.0 pass through at least 
four drift layers lying in a conic section outside of the 
central calorimeter. In both cases, muons are identi- 
fied as isolated tracks in the COT that extrapolate to 
track segments in one of the four-layer stacks. Two 
Cerenkov Luminosity Counters (CLC) in the forward 
region measures the average number of interactions per 
bunch crossing. 

The CDF trigger system is structured in three levels. 

Level 1 trigger uses informations from the COT, the 
calorimeters and the muon chambers. The eXtremely 
Fast Tracker (XFT) processor reconstructs the charged 
tracks using the hits from the axial layers of the COT 
with reduced resolution. Electrons, photons and jet can- 
didates are identified imposing the presence of energy in 
the single towers of the calorimeters above the thresh- 
old values. The value of the sum of the energy released 
on all the towers is used for the selections based on the 
total transverse energy and the missing transverse en- 
ergy. The transverse projection of the tower energies 
(see Sect.© are calculated with the assumption that the 
event primary vertex is located at z = 0. The missing 
energy at LI has poor resolution, due to a limited avail- 
able information and the need to make a fast decision, 
and is usually underestimated. 

At Level 2, the XFT tracks with > 4 hits in SVXII 
and pr > 2 GeV/c are reconstructed taking advantage 
of the additional information from the silicon detector. 

Algorithms can be applied to trigger electrons, pho- 
tons or jets: the energy for clusters of adjacent towers, 
and the information from the detectors of maximum ex- 
pansion of the shower in the electromagnetic calorime- 
ters are now available. 

Level 3 uses complete informations supplied from the 
various detectors. 

In particular, the tracking is completed executing the 
three-dimensional reconstruction of the trajectories in 
the volume \rj\ < 2, and more detailed algorithms recon- 
struct the energy in the calorimeters. The information 
provided by the cluster finding algorithm at the trigger 
level can be considered as a first-order jet reconstruc- 
tion. At this level, a looser time constraint enables to 
exploit the full detector segmentation for a better jet en- 
ergy and direction determination. 



3. Evaluation of the invisible particles' transverse 
momentum 

Since the discovery of the W boson, the missing 
transverse energy has been a crucial tool in the search 
for new phenomena at hadron colliders. Particle col- 
liding at hadron colliders have equal and opposite mo- 
menta; therefore, the total vector momentum sum in an 
event should be zero. The hard collision happens be- 
tween the partons of the proton (antiproton), where they 
can carry any fraction of the parent proton or antipro- 
ton momentum. However, in the plane transverse to the 
beam, the initial partons' momentum - and thus final 
ones as well - can be treated as zero to a good approxi- 
mation. Any transverse energy imbalance in the detec- 
tor may indicate that a particle left the detector without 
interacting with its material. 

Missing transverse energy, t^, is defined as the mag- 
nitude of the vector 



■ Z E t^ 
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where E l T are the magnitudes of transverse energy 
measured in each calorimetric tower i, and Hi is the unit 
vector from the interaction vertex to the tower in the 
transverse (x, y) plane. In order to maintain good reso- 
lution on both magnitude and direction of the missing 
transverse energy, the detector calorimeter must be de- 
signed to be almost completely hermetic, within the me- 
chanical allowance. 

The ¥t resolution is parametrized in terms of the 
total scalar energy deposited in the calorimeter (ZE T ), 
and is measured with minimum bias events, i.e. events 
collected requiring hits in both CLC counters, with no 
calorimeter requirements. No significant P r is expected 
in minimum bias events. A fit to these data yelds: 



o-(¥ Tx ) = -0.60 + 0.74 V^r 



(2) 



Similarly to the case of the missing transverse energy, 
it is possible to define a missing transverse momentum 
0t using the charged particle spectrometer, as the nega- 
tive vector sum of the charged particles momenta pr : 



= - Z Pt 

tracks 



(3) 



In events where only charged particles and unde- 
tected particles are produced, the 0t is highly correlated 
in module and direction to the undetected particle(s) 
momentum, and thus provide a way to measure their 
energy with potentially better resolution than £ T . The 
presence of quarks/gluons in the final state complicates 
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the picture. In the parton shower and hadronization pro- 
cess forming a jet, most particles produced are pions, 
with a smaller contribution of kaons. Due to isospin 
symmetry, roughly 2/3 of the energy of a jet will be 
carried by charged particles, which will be measured 
with both calorimeter and tracking chamber. The $t 
underestimates the undetected particle's energy because 
it does not take into account the energy carried by the 
neutral components of the jets, and for the same rea- 
son has a worse angular resolution than In addition, 
while calorimeter coverage is usually almost complete, 
the spectrometer coverage in collider detector is gen- 
erally far more limited. For these reasons, $j cannot 
substitute the role of the Mt as a tool for the measure- 
ment of momentum and direction of undetected parti- 
cles; still, the 0t can provide informations complemen- 
tary to those given by Er- 
in order to reconstruct the composite 0t observable, 
we use as inputs to its computation the charged parti- 
cle momentum reconstructed in the spectrometer. The 
basic track quality criteria for the calculation of 0t are 
the same used for the reconstruction of the primary ver- 
tex at CDF. Only tracks with 0.5 GeV/c < p T < 200 
GeV/c, \rj\ < 1.5 and \Z vtx \ <2 cm are used, where Z V f X is 
the closest approach distance of the track from the pri- 
mary vertex along the z axis. We then classify the tracks 
in four different categories, on the basis of the number 
of axial (N a c x 0T , N a s x vx ) and stereo (#£ or , N* vx ) COT 
and SVX layers that have at least 5 hits, together the 
X 2 of the track fit. The first category fulfills strong re- 
quirements on the track reconstruction in the COT; the 
second, third and fourth have decreasing requirements 
on COT compensated by requirements on the^ 2 and on 
the SVX hit layers. If the first category requirements 
fail, the second category is checked, then the third and 
finally the fourth. 

3.1. Data and simulated samples used 

We use for the following studies 5.7 fb" 1 of data, 
collected requiring the presence of large and the 
presence of two or more energetic hadronic jets. At 
Level 1 , the trigger requires at least one trigger tower 
with E T > 10 GeV, and a ¥ T > 28 GeV. At level 2, 
at least 2 calorimetric clusters with transverse energy 
above 3 GeV and \rj\ < 3.6 are required. The is 
recomputed with the additional available informations; 
again, $ T must be greater than 28 GeV. Finally, the 
Level 3 performs the complete ¥t determination, re- 
quiring it to be greater than 30 GeV. 

Events containing large ¥t can originate from non- 
collisional sources, such as cosmic or beam-halo muons 
passing the detector or noisy/dead calorimeter cells 



causing energy imbalance. These types of events are re- 
moved by requiring that the event observables indicate a 
presence of inelastic collision with large energy transfer, 
such as the presence of at least one high quality primary 
vertex in the collision. Additional requirements are also 
imposed to remove events consistent with beam-halo 
muons traversing the detector or those caused by noisy 
calorimeter cells. We require the corrected to be 
larger than 50 GeV, and we veto events containing at 
least one isolated electron or muon. We require the 
transverse momentum of the jet pr(i) (where i runs over 
the total number of jets) to be greater than 40 GeV for 
dijet events, and greater than 30, 20 and 15 GeV respec- 
tively for (i =1, 2), (i = 3, 4, 5) and (i > 5) for all other 
jet multiplicities. All jets must satisfy \rj\ < 2.4. With 
these selections, the QCD multijet production accounts 
for more than 95% of the data 11411 ; for this reason, we 
will mostly use the data themselves to show the feature 
of the QCD background. 

On the other hand, in order to describe the quali- 
tative features of physics processes giving rise to un- 
detectable particles in the final state, we use detailed 
Monte Carlo (MC) simulation interfaced to the detec- 
tor response simulation. For simplicity, we restrict our- 
selves to the modeling of the W/Z boson plus jets pro- 
duction, and top quark pair production. As an example 
of interesting signal, we use here the associated produc- 
tion of a Higgs (H) boson together with a W or a Z bo- 
son production WH/ZH, where the Higgs boson decays 
to bb and the Z (W) boson decays to vv (£v). The associ- 
ated production of a W or Z boson and jets is simulated 
using the alpgen Il5ll program, interfaced with parton- 
shower model from pythia 11611 . A matching scheme 
is applied to avoid double counting of partonic event 
configurations 1 17]. The samples are normalized to the 
inclusive cross sections 11811 . scaled by 1.3 to account 
for next- to-leading-order corrections. We model tt us- 
ing pythia with top quark mass equal to 172.5 GeV/c 2 , 
which is consistent with the current world best estimate 
of this parameter Jl9l |20|] and normalizing its contribu- 
tion to the tt NLO cross section 1 21]). Associated pro- 
duction of a Higgs (H) boson together with a W or a Z 
boson production WH/ZH is modeled with the pythia 
MC generator, assuming a Higgs boson mass equal to 
115GeV/c 2 . In order to test MC prediction for QCD 
events, we generate dijet and trijet events with pythia. 

For all the simulated samples above, the detector re- 
sponse is modeled by a geant 122] -based simulation of 
the CDF detector. 
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3.2. ¥j and $t distributions in QCD multijet events 



The generic production of strong-interacting parti- 
cles by quantum chromodynamics interactions between 
quarks and gluons is the most probable source of 
hadronic jets in the particle colliders. Because of this, it 
is also the greatest source of backgrounds for searches 
which do not require leptons in the final state. QCD 
multijet production does not include undetectable parti- 
cles, apart from small contributions from heavy quarks 
semileptonic decays, which produce neutrinos. Missing 
transverse energy in these events, if correctly measured, 
is thus expected to be zero. This allows to significantly 
reduce the QCD multijet background for signatures in- 
cluding undetectable particles. Nevertheless, calorime- 
ter resolution as well as non linearity and non hermetic- 
ity have to be taken into account, resulting in a smearing 
of the ¥t distribution. As a result of this, the largest part 
of the events which pass a P r > 50 GeV requirement, 
which is loose enough to be used with low final momen- 
tum signatures (see Sect. [B, is mainly constituted by 
QCD multijet. Furthermore, this contribution is difficult 
to model, especially in topologies consisting of large 
jet multiplicities. For this reason, analysis looking for 
^r+jets generally take advantage of data-driven tech- 
niques to model QCD contribution. In Fig. [T]are shown 
the ¥t distributions for events with ¥* T > 50 GeV, for jet 
multiplicities ranging from 2 to > 7; as can be seen, the 
QCD multijet is concentrated to lower values, while 
MC simulations of tt, W+jets and Z+jets have, except 
in full hadronic decay modes, real P r due to neutrinos. 

Again, the data themselves are taken as a represen- 
tation of the QCD background. The distribution shows 
that physics processes giving rise to high-p T neutrinos 
display a P r distribution with a pronounced tail even 
after a given cut on the magnitude of the missing trans- 
verse energy - as expected. 

Figure[2] shows the $t distribution for the same 
physics processes, after the same event selections. It 
can be noted that the 4t distribution for events contain- 
ing real neutrinos in the final state tend to peak close to 
the IV values for the same processes. On the other hand, 
the data - representing QCD multijet production - peaks 
at much lower values, in a manner compatible with the 
interpretation that $t in QCD events arises from ran- 
dom fluctuations of the charged particle component of 
each jet in the event, and that they tend to cancel out in 
the simple dijet topology. The bottom line is that the 
$t clearly shows a notable separation power between 
events where ¥t arises because of resolution effects, 
and physics processes where is expected to signal 
the passage of undetectable particles. 
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Figure 1: The distribution of ¥t for data and major Standard Model 
processes giving rise to neutrinos. The events are selected with ¥t > 
50GeV and number of jet Nj et varying from 2 (top-left plot) to 7 or 
more (bottom-right plot). All distributions are normalized to unit area. 
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4. Angular correlations between ^and 0t 
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Figure 2: The distribution of $t for data and major Standard Model 
processes giving rise to neutrinos. The events are selected with > 
50GeV and number of jet Nj et varying from 2 (top-left plot) to 7 or 
more (bottom-right plot). All distributions are normalized to unit area. 



For events with real undetected particles, £j is a 
good measurement of the vectorial sum of those par- 
ticles four-momenta , and so is the which thus will 
tend to be parallel to the £j- However, in QCD events 
the two variables have very different origin, which is 
reflected in their different expected angular correlation. 
We will describe this aspect in details, starting from the 
simplest topology. 

4.1. Dijet events 

For QCD events, the energy conservation requires 
that the vector sum of the jet transverse energies 
amounts to zero. In dijet events, the 2 jets will have 
the same magnitude of transverse energy, and will come 
out back- to-back in azimuthal space. The mismeasure- 
ment of jet energies makes the Mj align to the jet with 
less measured energy. The 0t could be present in some 
amount too, but for different reasons. In the jet fragmen- 
tation and hadronization processes a certain number of 
charged particles are produced inside each jet, and they 
will be detected by the tracker. The intrinsic fluctua- 
tions inherent to the parton shower process will result 
in large fluctuations in the fraction of energy carried by 
charged particles inside a jet. The fraction of charged 
particles for a jet is completely independent from the 
measured jet energy in the calorimeter. Since is more 
likely for a jet's energy to be undermeasured than to 
be overmeasured, the Mj will always be aligned to the 
under-measured jet, while the fj will mainly point to 
the direction of jet with less energy being carried by re- 
constructed charged particles. As a net result, in QCD 
dijet events the Mj and 0j directions will be mainly cor- 
related or anticorrelated. A schematic representation of 
the argument above is given in Fig. [3] 

This particular topology is reflected in the distribu- 
tion of the azimuthal distance between the two vectors, 
A(p(£ T ,0 T )- As shown in Fig. |4l events containing un- 
detected particles like neutrinos are concentrated near 
0, while QCD dijet events have almost equal chance to 
populate the region around or around n. This spe- 
cial distribution allows to effectively suppress the QCD 
contribution in samples with ¥t and jets. Dijet (and 3 
jets) events are important because have the most sim- 
ple kinematic, so that they can be also studied through 
QCD MC. Fig. [5] shows the A0(i? r ,^ r ) distributions 
for data and simulated samples of QCD and events with 
neutrinos. The presence of b quarks in the final state 
is expected not to alter significantly the situation, be- 
cause their semileptonic decays would give rise to low- 
energetic neutrinos. 
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Figure 3: Schematical representation of the energy and momentum 
flow in a QCD dijet event. In both plots, Ej is due to calorimetric 
mismeasurement. The $t, which is due to statistical fluctuation in the 
vectorial sum of charged particles' momenta in each jet, can arise with 
equal probability in the correctly-measured jet (top) or in the under- 
measured jet (bottom). 
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Figure 4: The distribution of A(P(Et,$t) data and major Standard 
Model processes giving rise to neutrinos. The events are selected with 
¥t > 50GeV and number of jet Nj et varying from 2 (top-left plot) to 
7 or more (bottom-right plot). All distributions are normalized to unit 
area. 
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Figure 5: The A0(^V,/r) distribution in events with 2 (left) and 3 
(right) jets. Lines represents QCD, Higgs+W/Z and MC, while dots 
represents data. All plots are normalized to unit area. 



The effects of fluctuations in the charged fractions of 
the jet energies can be quantified defining the charged 
fraction as follows: 



chgd _ 1 ^Jet trac ki r i i\ ... 
Pt ~ jet ^ ' 

P T 

In QCD dijet events the requirement of large ¥t 
means that the energy of one of the two jets energy will 
be largely mismeasured. After ordering the jets accord- 
ing to their Ej starting from the one with the largest en- 
ergy, the Mt will be aligned to the second jet. In events 
where Mt and $t are correlated we expect the p C j gd of 
the 2 nd jet to be sistematically smaller than in events 
where they are anti-correlated. The opposite effect will 
be present for the 1 st jet. This can be observed by look- 
ing at the distribution of p C j gd of the I st and the 2 nd jet, 
separately for events where Mt and 0t are correlated or 
anti-correlated. As can be seen in Fig. [6] the expected 
results are confirmed, which supports the overall inter- 
pretation of the mechanism behind the formation of the 
A(p(M T , 0t) distribution for QCD events. 



4.2. Large jet multiplicities 

In energetic QCD multijet events, large part of the 
energy is carried by the two most energetic jets. These 
two jets will tend to be roughly back- to-back, and both 
Mt and 0t will have high probability to lie on their axis, 
while additional jets will arise because of radiation from 
the initial/final state particles. For these reasons, we ex- 
pect the exposed A(/)(Mt, $t) distribution features to be 
recognizable also for QCD multijet events. This is con- 
firmed from the analysis of both data and MC with ex- 
actly 3 jets, as is shown in Fig. |4]and[5] The logical step 
can be extended to higher multiplicities. In fact, for a 
typical QCD multijet event with 4 or more jets, more 
than 60% of the total energy is carried by the two most 
energetic jets. The distributions of data shown in Fig. [4] 
confirm the validity of this approach up to events with 
7 jets or more, that is where CDF statistics run out. The 
use of MC to cross-check, as for 2 and 3 jets events, is 
not feasible for higher multiplicities; in facts, because 
of the high production rate for QCD at a hadron col- 
lider and the large statistics needed in order to describe 
these processes adequately in an analysis of several in- 
verse femtobarns of data, the simulation of an accept- 
able amount of QCD events is prohibitive. 

5. Study of ¥<t and jf T as a function of the collider 
instantaneous luminosity 

The varying instantaneous luminosity at hadron col- 
liders poses a unique challenge to the online selection 
(triggering) of events characterized by large missing 
transverse energy and jets in the final state. In fact, 
larger instantaneous luminosity means more collisions 
per bunch-crossing. The most common collisions will 
give rise to quarks or gluons in the final state, and thus 
jets in the detector. These jets can add up to the miss- 
ing transverse energy from the hard scatter, or can pro- 
vide additional missing transverse energy if the energy 
of any or more than one of those is mismeasured in the 
detector. 

In both cases, the pile-up of multiple scatters will 
constitute an additional background to the missing 
transverse energy plus jets signature. Typical solutions 
proposed to cope with the unsustainable large trigger 
rates at increasing luminosity, have been to tighten the 
requirement of the absolute value of the missing trans- 
verse energy, or to tighten the requirement in the jet 
transverse energy, or both. While allowing to control 
the trigger rates, this procedure often reduces drastically 
the yields of the interesting signal. 
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Figure 6: The pf gd distributions for the 1 st (left) and 2 nd (right) jet in QCD MC, while A<f>(M Tl p T ) is in (0, tt/2) or (tt/2, tt); the dots shows the 
same distributions for data. MC expectations are normalized to data. 



We argue in this section that the spectrometer can be 
helpful in suppressing the multiple scatter pile-up back- 
ground in a similar fashion as proposed in suppressing 
the single- scatter QCD background. 

To study the dependence of the pfj on the instanta- 
nous luminosity, we choose as a proxy the number of 
primary collision vertices N vtx as reconstructed in the 
spectrometer. While this constitutes an approximation, 
since for any vertex reconstruction algorithm efficiency 
drops progressively as the luminosity increases, it is suf- 
ficient to understand the qualitative features of luminos- 
ity dependence. 

Analyzing the data for different numbers of interac- 
tion vertices is possible to appreciate the pile-up rejec- 
tion power of a cut on $ T . This can be particularly in- 
teresting in a LHC context, where the far higher instan- 
taneos luminosity makes extremely important to control 
this source of background. Fig. [7] shows the IV, $t and 
A<P($t,$t) distributions dependence from the number 
of vertices for the whole data sample (/. e. for any num- 
ber of jets). As can be seen, after fixing a certain online 
cut on IV, the offline ¥t has only a mild dependence on 
the number of primary vertices. On the other hand, the 
pile-up events tend clearly to populate the region with 
QCD-like characteristics, i.e. at low-P^r and simmetri- 
cally in the A(p(£ T , T ) distribution. 



6. A A0(^t> /jO-based QCD data-driven model 

We stressed that the peculiar way QCD populates the 
A(p(£ T ,0 T ) distribution is by itself a powerful instru- 
ment to reduce the QCD background by a factor of two 
when cut at 7r/2, with minimal signal loss. Its imple- 
mentation in more complex methods, like multivariate 
analysis, would maximize its impact. 

We argued, assisted by a MC simulation, that in the 
simplest case QCD populates the region where $t is cor- 
related or anti-correlated with $ T depending on the frac- 
tion of energy of the jets carried by charged particles. 
The fluctuations in the vectorial sum of charged parti- 
cles' momenta in a jet, and in the amount of transverse 
momentum carried by those, derives from fluctuations 
in the parton shower and hadronization non-perturbative 
process. On the other hand, the ¥t in QCD events ap- 
pears preminently as a consequence of the fluctuation 
in the jet energy measurement in the calorimeter, or in 
the peculiar topology where one or more of the jets in 
the event falls in an uninstrumented region of the de- 
tector. The reasons for the appearance of the IV and of 
the $t in QCD events are completely uncorrected. An 
angular correlation is on the other hand introduced by 
the fact that in the simple event topology of QCD di- 
jet events, ¥t and $t have to be aligned with either one 
jet or the other. As demonstrated in a previous section, 
this topology is only mildly altered in the circumstance 
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Figure 7: From top to bottom: the Fr, $t and A^(Et,^t) distribu- 
tions in events with varying number of interaction vertices. It is clear 
how the latter two plots show far better discrimination power toward 
the isolation of pile-up events with respect to the traditional cut on 
missing transverse energy, and can thus be used to control more effi- 
ciently the pile-up background. All plots are normalized to unit area. 



where QCD events appear with larger jet multiplicities 
due to radiation. 

An important consequence of the peculiar angular 
correlation between $ T and $t is that QCD is supposed 
to populate symmetrically the A^{£ t ,0t) distribution, 
with the symmetry centered at n/2. By rejecting the 
events with A(p(£ T ,$ T ) > n/2 and defining as signal 
region the one with A<p(Er,pr) < 7r/2, one can infer 
the contribution of QCD in the signal region in a data- 
driven fashion by counting the events populating the re- 
jected region, i.e. A(p(£ T ,0 T ) > n/2. The robustness of 
this method has been tested by varying the jet multplic- 
ity, by placing cuts on the angular correlation among $ T 
and the jets direction, varying cuts, etc. 112311 . In all 
instances, by using MC simulation the above assump- 
tion has been found true to a level of ~ 20%. The exact 
value should be determined by properly defining a "con- 
trol region" for the analysis. 

Following the same reasoning, it is legitimate to as- 
sume that QCD events in the A(f)(M T ,0 T ) > n/2 could 
be used to model the kinematic features of QCD events 
in the complementary region. The systematic uncertain- 
ties associated with the simulation of high multiplicity 
QCD jet production are very large. For these reasons 
a reliable data-driven modeling method for this cate- 
gory of events is of fundamental importance in collider 
physics analysis. This has been done in J2|,|24j] for dijet 
signatures and in HHSl with multiplicities up to 10. In 
all instances, a general good agreement has been found 
for a large number of the kinematical and topological 
distributions. 

7. Outlook on the usage of jf T at the ATLAS and 
CMS detectors 

The ATLAS and CMS calorimeter performance are 
comparable or mildly superior to the one of the CDF 
detector used in this study. This means that QCD multi- 
jet production is still a serious obstacle to physics mea- 
surement and searches in the P r + jets signature. On 
the other hand, the coverage of the ATLAS and CMS 
charged particle spectrometer is far superior to the CDF 
one. Both arguments lead the authors to believe that the 
application of the techniques presented here to these ex- 
periments would be straightforward, and bring notable 
advantages. In particular, the authors call for a revisita- 
tion of the current ATLAS and CMS strategies of tight- 
ening the ¥t requirements to cope with the increasing 
collider instantaneous luminosity, and suggest and alter- 
native way to control the QCD background. The authors 
are aware that CMS uses already the spectrometer in a 
"particle flow" reconstruction in order to increase the P r 
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and jet energy resolution; still, improving the resolution 
on those observables reduces the QCD contribution to 
IV + jets samples but does not fully eliminate it. We 
argue that even after the use of particle flow techniques 
the $t and its correlation with ¥t will help to suppress 
QCD and to model its contribution. It is foreseeable that 
future upgrades of the ATLAS and CMS trigger systems 
could allow the computation of a $t "primitive" online, 
in order to contain the and P r + jets trigger rates 
to reasonable levels, without having to raise P r cut to 
the point where the loss of interesting signals become 
unacceptable. 

8. Conclusions 

We introduced here a new observable that can be used 
to infer informations on the four-momentum of invisible 
particles passing trough collider detectors. The paper 
focuses on collisions where invisible particles appear 
together with quarks or gluons in the final state, where 
improvements in the identification techniques of invisi- 
ble particles is crucial in order to suppress the otherwise 
dominant QCD background. 

As the technique uses only the spectrometer, as op- 
posed to the traditional use of the calorimeter, the instru- 
mental resolutions are completely uncorrected. In par- 
ticular, the resolution on the $t measurement depends 
almost solely on the inherent fluctuation on the num- 
ber and transverse momentum carried by the charged 
hadrons of the jets in the event. 

We have further shown the importance of the evalua- 
tion of the angular correlation between T and Mj, as a 
valid tool to reduce the contribution of the QCD multijet 
background in t^+jets signatures, studying the proper- 
ties of ^--related observables in collider data and in MC 
simulations. By taking advantage of the peculiar prop- 
erties of $t, we developed a data-driven way to under- 
stand the QCD contribution to the ¥t + jets sample with 
relatively large precision, and suggest that in a similar 
fashion one can reproduce with similar level of detail 
the kinematic properties of QCD multijet events. The 
techniques developed here have been already success- 
fully used in a broad range of a naly ses at the Tevatron 
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